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Abstract 

An operator deformed quantum algebra is discovered exploiting the quantum Yang-Baxter equa- 
tion with trigonometric R-matrix. This novel Hopf algebra along with its g — + 1 limit appear to 
be the most general Yang-Baxter algebra underlying quantum integrable systems. Three different 
directions of application of this algebra in integrable systems depending on different sets of values 
of deforming operators are identified. Fixed values on the whole lattice yield subalgebras linked to 
standard quantum integrable models, while the associated Lax operators generate and classify them 
in an unified way. Variable values construct a new series of quantum integrable inhomogeneous 
models. Fixed but different values at different lattice sites can produce a novel class of integrable 
hybrid models including integrable matter-radiation models and quantum field models with defects, 
in particular, a new quantum integrable sine-Gordon model with defect. 

PACS numbers 02.30.Ik, 02.20.Uw, 05.50-Fq, 03.65.Fd 

Key Words: Operator deformed quantum algebra; unifying scheme for quantum integrable systems; 
new quantum integrable models: inhomogeneous, matter-radiation, sine-Gordon with defect. 

1 Introduction 

Quantum integrable systems (QIS) are nonlinear, interacting and nonperturbative quantum many 
body or field models having large symmetries with mutually commuting set of conserved operators 
{Cn}, n = 1, 2, . . .. Such models usually allow exact solution of the eigenvalue problem through the 
algebraic Bethe ansatz method, for all conserved operators including the Hamiltonian of the system. 

Every lattice QIS is represented by its own Lax operator Lj{\) at each lattice site j = 1,2, 
and for the integrability it should satisfy quantum Yang-Baxter equation (QYBE) 

R{X - iJi)Lj{\) (g) Ljifi) = {Lj{fi) I){Lj{X) ® I)R{X - n), j = 1, 2, N (1.1) 

together with the ultralocality condition Lj{X) Lk[fi) = (L^dj,) I){Lj{X) I), at / j. These 
two relations combined together lead to several important consequences. 
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i) They provide sufficient condition for quantum integrability. Defining monodromy matrix T{X) = 
Ylj Lj{\) on the whole lattice, one shows that the same QYBE holds also for this global object : 

R{\ - /x)r(A) Till) = (Till) (g) /)(r(A) ® I)R{X - /x). (1.2) 

Taking trace from both sides one derives the trace identity: [r(A),r(/u)] = 0, for r(A) = Tr{T{X)), 
yielding the crucial integrability condition [Cn,Cm] = 0, through t(A) = S„ C„A" 

ii) QYBE (1.2) plays important role in exact solution through the algebraic Bethe ansatz method by 
providing commutation relations between diagonal (generator of conserved operators) and off-diagonal 
('creation' / 'annihilation' operators ) elements of r(A)-matrix. 

iii) Finally, QYBE (1.1) defines the commutation relations between operator-elements of the quan- 
tum Lax operator (Yang-Baxter (YB) algebra), specified by the quantum i?-matrix, a c- number matrix 
which fixes the the structure constants of the algebra. Interestingly, the very fact that the same QYBE 
is valid also at the global level (1.2), tells that this YB algebra must be a Hopf algebra A with defining 
properties: 

coproduct: A : A ^ A<Si A , product: m: A<Si A ^ A, unit: e : k ^ A, counit : e : A ^ k 
and antipode (inverse): S: ^ — ^ ^. 

For some integrable systems when the last property is absent the algebra becomes simply a bi- 
algebra (algebra -|-coalgebra). The underlying YB algebra classifies the quantum integrable models 
into three major classes : rational, trigonometric and elliptic, distinguished by such solutions of the 
associated quantum i?-matrix. Here we will focus on the first two classes. 

2 Known Yang-Baxter algebra 

The simplest YB algebra is given by the standard Lie algebra sl{2) : [S+ , S~] = 2S^, [S^, S^] = ±S^. 
This is a cocommutative Hopf algebra generated by the quantum i?-matrix: 

-^11 = -^22 = fl(-^)) -^12 = -^21 = K^)j -^21 = -^12 = c, (2.3) 

having solutions i?*""* in rational functions of the spectral parameter: a(A) = X + a, 6(A) = A, c = a. 
This algebra is linked to the well known quantum integrable models [1] like 1) Nonlinear Schrodinger 
(NLS) equation and the lattice NLS, 2) xxx-spin-i chain 3) rational Gaudin model etc. 

A more general (and more interesting) YB algebra is associated with the trigonometric 
matrix solution of (2.3) as a(A) = sin(A-l-a), 6(A) = sin A, c = sin a. and is given by the celebrated 
quantum algebra slq{2) [Sq,S~] = ^'"^^."^^ , [S'^,S'^] = ±5"^, where q = e*" is the deformation 
parameter. FRT construction [2, 3] gives an elegant way of obtaining this quantum algebra together 
with its non-cocommuting Hopf algebra properties and the explicit form of the associated (2 x 2) Lax 
operator. slq{2) was shown [3] to be the underlying YB algebra of well known quantum integrable 
systems like 1) sine-Gordon (SG) and lattice SG model, 2) xxz-spin chain etc. 

We however observe that there exists a series of other important quantum integrable models 
of diverse nature having completely different Lax operators but associated with the same quantum 
matrix: i?*''*^ or BJ"""^. Few such examples are the quantum Toda chain belonging to the rational class, 
and the relativistic Toda chain, Liouville model, derivative NLS etc. belonging to the trigonometric 
class. Therefore in spite of the far-reaching success of the quantum algebra 5/^(2) it is pertinent to 
ask 
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i) Does it exhaust all possible YB algebras associated with the trigonometric i?-matrix ? ii) What 
are the YB algebras underlying other integrable models ? iii) Can they all be obtained from 5/^(2) 
or sl{2) ? iv) Can there be any deep relation among diverse integrable models sharing the same 
quantum i?-matrix, despite of having completely different Lax operators ? 

In search for answers to these fundamental questions, we conjecture that there must exist a more 
general YB algebra and the related quantum Lax operator associated with the same i?*''*^-matrix 
(and its g — >^ 1 limit linked with i?''"*-matrix) , which should cover all quantum integrable models with 
(2 X 2) Lax matrix. It should include the known quantum algebras (and its g — >^ 1 limit) along with 
related models as particular cases, while having the freedom to generate YB algebras underlying all 
other known and new quantum integrable models. 



3 Operator-deformed quantum algebra 

Supporting our conjecture we discover a conceptually novel YB algebra, an algebra deformed not only 
by usual parameter g, but also by some operators. Our idea is to follow closely the FRT construction 
for slq{2) and the associated Lax operator [2], but deform the related Borel subalgebras further by 
operators c^, a = 1, 2, resulting a new Lax operator 

^anclU 1^ 2sina4 ^c+e-^°^'+r'c2-e^"^' ' ^ ' ^ ' 

while the i?*"^(A)-matrix remains the same. Inserting (3.4) together with in QYBE (LI) we derive the 
novel operator-deformed quantum algebra, we call Ancestor algebra [4], defined through the relations 

S-] = (m+ sin(2as^) - iM' cos{2as^)) -J—, [s^, s^] = ±s^, [M^, •] = (3.5) 

where the deforming operators = ^{cfc2 ± cJ~C2') and all operators c~±, a = 1,2 are mutually 
commuting as well as central (commuting with all other generators of algebra (3.5)). We believe that 
(3.4) is the most general form of the Lax operator, which can generate 2x2 Lax operators of all 
integrable systems belonging to this class and the operator-deformed quantum algebra (3.5) is the 
most general YB algebra allowed by the simplest i?*^*^-matrix. This ancestor algebra exhibits the 
following unique and distinguishing properties. 

1). The operator deformed quantum algebra (3.5) is a bialgebra (algebra + coalgebra) with the 
Hopf algebra properties: 

i) Coproduct A : 



^(s^) = s+ Ci g + C2 g"' s+, A{Sg ) = Sg (g) c^q + cfq'" ® Sg , 
A{s^) = / (g) + (g) J, A(c^) = (g) c^, a = 1,2 

Note that unlike deforming parameter q, deforming operators and hence have nontrivial 
coproduct. 

ii) Antipode S: 

S{s-) = -(c+)-ie--^%-e-^'(c+)-\ S{s^) = -{c^r'e-"''' S+e"''\c^r\ 
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In a QIS if any c^=0, the antipode disappears turning the Hopf algebra into a bialgebra. 
iii) Counit e: 

6(.±) = 0, 6(e±-^') = l, 6(C±) = C± 

For product m one can take the formal definition of multiplication in the algebra, while unit rj may be 
defined through the unital element 1 as r/(,^) ^1. 

2). Note that algebra (3.5) is deformed by parameter q as in the usual quantum algebra, but 
additionally it is deformed also by central operators Af^, prompting our definition of this deformed 
algebra as a novel operator deformed quantum algebra. The deforming operators being central, they 
act simply as multiplicative parameters, which however can take positive/negative or zero values, 
generating different possible subalgebras as described below. It should be noted that for quadratic 
algebra (3.5) more sophisticated construction like quantum double, universal R-matrix or even finite- 
dimensional representations might be problematic and could be an interesting future problem to 
explore. However this does not create any problem in applying it as a YB algebra for generating and 
classifying QIS, which is the main purpose of the present investigation. Generically this algebra allows 
infinite-dimensional representation, which may be given through canonical operators [u, p] = i as 

s3 = 5+ = e-'Pg{u), s- = g{u)e'P (3.6) 

where ^ 

g'^{u) = ^ . 2 (« + sina(s - ■u)(M+ sinQ;(n -|- s -|- 1) - iM~ cos a{u + s + 1))) (3.7) 

with K and spin variable s being arbitrary constants. We drop hat from the deforming central operators 
which act only multiplicatively and can also have zero values. 

The integrable ancestor model, as obtained using (3.6, 3.7) in its Lax operator (3.4), represents 
a generalized lattice sinc-Gordon (LSG) model. It reduces to the known LSG [5] in the absence of 
deforming operators (i.e. at = 1, M~ = 0), when we also recover from (3.5) the well known slq{2) 
and from (3.4) the known Lax operator of [2]. In general (3.5) can yield nine major subalgebras for the 
range of values > 0, < 0, = (up to scaling), for each of the sectors M" = 0, >, < . Among 
these subalgebras M+ > 0, M~ = yields the well known slq{2), while M+ < 0, M~ = sector 
gives the corresponding noncompact case. All other cases give a variety of subalgebras describing 
underlying YB algebras of important quantum integrable systems, including new ones. 



3.1 Operator-deformed algebra at q ^ 1-limit 

Consider now q ^ 1 limit of the above operator-deformed algebra, when the related structures reduces 
to their corresponding rational limit: Quantum i?*''*^-matrix goes to the rational one i?*""*. The 
operators s^,b = 1,2,3, — c^'^,a =1,2 (since they are also q-deformed!) and the spectral 

parameter ^ ^ X, while the ancestor Lax operator (3.4) reduces to its rational form 

Trat/\\_( c\{X + s'^) + c\ S~ \ 

Finally algebra (3.5), with deforming operators — m^, reduces to another novel operator- 
deformed (but g — > 1) bi-algebra 

[s+,s-]=2m+s^ + m-, [s^,s^]=±s^, [m^,-]=0 (3.9) 
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with central operators m"^ = c^c^, rh~ = c}c2 + C1C2, as expressed through mutually commuting 
central operators cf, a = 0, l;i = 1,2. Interestingly it still exhibits non-cocommutative coproduct 
given by A(s'*") = s"*" + s"*" ® 62, A(s~) = C2 (8) s~ + Cj*, etc. Note that for nontrivial rh~ 
this is a new operator-deformed Hopf algebra with nine distinct subalgebras with > 0, < 0, = 0, 
for each of the combinations from = 0, > 0, < 0. The case > 0, m~ = yields the standard 
Lie algebra sl{2), while m+ < 0, rn^ = yields its noncompact variant sl{l, 1). The rest of the cases 
gives the freedom of describing other integrable models. 

3.2 Application to quantum integrable systems 

Since only the ancestor Lax operator (3.4) contains the deforming operators but not the associated 
i?*"^-matrix, we can construct Lax operators of all integrable systems belonging to the same trigono- 
metric class together with their YB algebras, by suitable choices of the deforming operators. Similarly 
at g ^ 1 limit we would get the corresponding result for the models belonging to the rational class 
all sharing the same ii''"*-matrix. We identify three types of integrable models that we can construct 
following our scheme, depending on different range of values of the deforming operators. 

i) For fixed values of the deforming operators on the whole lattice we get from (3.5) its different 
subalgebras which interestingly represent the underlying YB algebras of original integrable models. 
At the same time the ancestor L^"|' (^)-operator (3.4) and its rational limit L'^^ric (3-8) generate in a 
unified way the representative Lax operators of known as well as new quantum integrable models as 
R(xxz) : L^anc i)xxz — spin chain , ii) sine- Gordon (lattice +field) 

, Mi)derivative NLS (lattice -|-field ) , TO)q-bosonic , 'u)Relativistic Toda chain 
■uz)Liouville (lattice -)-field) 'yM)massive Thirring (bosonic) 

And at g ^ 1 (rational limit) 
R{xxx) : L™*^ =^ i)xxx — spin chain , ii) NLS (lattice +field) 
in) simple lattice NLS iv) Toda chain 

This systematic generation of quantum integrable models from a single ancestor model also answers 
to some intriguing questions raised above. Our finding shows clearly that all these integrable models 
are related deeply to each other as descendants of the same ancestor and therefore they share the same 
inherited R-matrix. Moreover their underlying YB algebras are obtained as different subalgebras of 
the same ancestor algebra. 

ii) For variable values of the deforming operators on the lattice however we obtain a new series of 
quantum integrable inhomogeneous models. 

iii) Choosing different but fixed values of the deforming operators at different sectors of the lattice, 
we can combine different interacting integrable models and generate a novel class of integrable hybrid 
models. 

We present below details of the YB algebras as subalgebras of our ancestor algebra and its major 
applications in quantum integrable systems as mentioned above. 

4 Generation of integrable models and underlying YB algebras 

Starting from the ancestor Lax operators (3.4) and (3.8) and the related operator deformed algebras 
(3.5) and (3.9) we generate quantum integrable models, known as well as new, together with their 
underlying YB algebras in a unified way. 
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4.1 q-deformed algebra and related integrable models 

We first focus on models belonging to the trigonometric class with i?*"^. 

1) At q-deformed but operator-undeformed limit: cf = C2 = = C2 = I giving M+ = 1, M~ = 0, 
our general algebra (3.5) reduces to the known quantum algebra slq{2) and the ancestor Lax operator 
(3.4) to the known L{X)- operator of [2]. Consequently the spin-^ representation would generate 
xxz-spin chain [3], while realization (3.7) reproduces g{u) = 2 sin a cosq;(2'u + l)]^ related to the 
known lattice SG model [5]. 

2) M+ = —1,M~ = case reduces (3.5) to the corresponding noncompact quantum algebra 
slq{l, 1) and reproduces the related integrable model, e.g. g-deformed Buck-Sukumar(BS) model. 

3) Taking the deforming operators as cf = (c^)~^ = p'^, we get p, g-parameter-deformed Hopf 
algebra glp^q{2), which is represented by the same slq{2) algebra, but is p"^ -deformed as a co-algebra. 
This YB algebra is linked with the quantum integrable Ablowitz-Ladik model. 

It is worth noting that, such an algebra is believed to be obtainable only through twisting trans- 
formation with an external parameter p = e'^ [6]. Under such a transformation the i?-matrix is also 
deformed to a twisted one: i?*''*^(A, ^)-matrix. However, contrary to this belief we obtain the same 
bi-algebra without changing the i?*"^(A)-matrix, since in our case the deforming operators replaces 
the role of twisting parameter as mentioned above. 

4) Choosing the values of the operators as M+ = sin a, M" = icosa, compatible with = C2 = 
1, c]" = —iq, C2 = ^ one gets the g-bosonic algebra 

['iPq,N]=^q, [,pl,N] = -^l, [V'^Vj] =cos(a(2Ar + l)), (4.10) 

by assuming s'^ = tpq, s~ = V'g; = N. This opens up an interesting possibility of constructing 
quantum integrable models involving g-boson [7]. In this case (3.7) simplifies to g^{u) = [—2u]q 

A g-boson can be mapped into a standard boson tp on the lattice with commutation relation 
= sis ijjq = ^(277^^^)^' ^ — V'V'' • This can construct from the above q-bosonic model an 
exact lattice version of the quantum integrable derivative NLS (with Lax operator obtained directly 
from (3.4)), which at the continuum limit A — > yields the quantum DNLS model [8] represented by 
the equation iijjt = ipxx — ■j('0^V')'0x- This QFT model in the N-particle sector is equivalent to the 
interacting Bose gas with derivative (5-function potential [9] . 

5) For another choice M+ = M~ = 1 we derive a novel exponentially deformed algebra 

which turns out to be the underlying YB algebra of the quantum integrable exact lattice Liou- 
ville model, with its Lax operator obtained from (3.4) with realization (3.7) reducing to g{u) = 

— ^^.^^ ' . At the continuum limit it gives the quantum Liouville field model given by the equation 

utt - Uxx = e"'' 

It is important to note that the present values of remain same even with c]" / 0, which 
would give the same algebra and the same realization, but a different Lax operator. This is indeed an 
intriguing possibility of constructing different useful Lax operators for the same model, in a systematic 
way. For example, a second, Liouville Lax operator we can construct here so easily recovers that 
invented by Faddeev through many innovative tricks [11]. 

6 ). In a similar way the particular case M=t = can be achieved with different sets of choices: 

i) c~ = 0, = 1, a = 1, 2 ii) = 0, = —cl = 1, Hi) cf = 1, rest of c's = 0, 
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all of which lead to the same algebra 

[s+,s-] = 0, [s^s^]=±s^. (4.12) 

However, they may generate different Lax operators from (3.4), corresponding even to different models, 
though with the same underlying algebra. In particular, case i) leads to the light-cone SG model, 
while ii) and iii) give two different Lax operators found in [12] and [13] for the same relativistic Toda 
chain. Since here we get g{u) =const. , interchanging u —ip,p —iu, in (3.7) yields simply 
= —ip, = ae^**, which generates a quantum integrable discrete-time or relativistic quantum 
Toda chain [12]. 

Remarkably, all above algebras are obtained as subalgebras from our ancestor algebra, but can not 
be obtained from slq{2), without invoking singular limits. Note also that all the descendant models 
listed above have the same trigonometric i?*'^'*^-matrix, inherited from the ancestor model and similar 
is true for the rational class, as wc will sec below. This unveils the mystery why a wide range of models 
found to share the same R-m&trix, while their L-operators are obtained from the same ancestor Lax 
operator (3.4) and their underlying YB algebras from the same ancestor algebra (3.5) at various 
reductions. 

4.2 q-undeformed algebra and related integrable models 

We focus now on the g — > 1 limit of algebra (3.5), which would still be operator deformed through 
rh^. It is interesting to find that the bosonic representation (3.7) at this limit reduces to a generalized 
Holstein-Primakov transformation (HPT) 

= s-N, s+ = go{N)i/;, s' = ^p^go{N), g^{N) = + m+(2s - N), N = V'V- (4-13) 

with the central deforming operators acting multiplicatively. This is an exact realization of (3.9), asso- 
ciated with the Lax operator (3.8). This rational ancestor model, representing a quantum integrable 
generalized lattice NLS model, can generate in a systematic way all other integrable models of the 
rational class through different choices for the values of the deforming operators. This constructs at 
the same time the Lax operators of the models as reductions of (3.8) and their underlying YB algebras 
as subalgebras of (3.9). 

1) When the deforming operators vanish with the choice m"*" = l,m~ = 0, (3.9) gives clearly the 
standard 5^(2) and for the spin ^ representation one recovers the xxx spin chain [3]. 

2) A bosonic realization of the same algebra simplifies (4.13) to the standard HPT and (3.8), 

reprodiicing the well known lattice NLS model [5]. 

3) . On the other hand for m"*" = —l,m~ = 0, (3.9) reduces to su{l,l) and gives the related 
integrable models, one of which is the integrable BS type model constructed in sect. 7. 

4) . For the complementary choice = 0,m~ = 1, (3.9) reduces to a non semi-simple algebra 

[s^,s']=Ts^, [s+,s-] = l (4.14) 

and (4.13) with go{N) = 1 gives a direct bosonic realization s'^ = ip,s~ = tjj^ , = s — N. Using this 
reduction in (3.8) we discover another quantum integrable simple lattice NLS model [14]. 

5) A trivial choice m"^ = reduces (3.9) to the same algebra (4.12) 

s~] = 0, [s^ s^] = ±s^ (4.15) 



7 



and therefore we can take the same reahzation = —ip, = ae^" found for the relativistic 
case, but now with ancestor Lax operator (3.8) and i?''"*-matrix. This gives quantum integrable 
nonrelativistic Toda chain [I]. 

Thus we have shown how the ancestor model (3.8) with (4.13) can generate quantum integrable 
models, all sharing the same rational i?-matrix in a unified way and their YB algebras like sl{2), 
(4.14), (4.15 ) are obtained as subalgebras from operator-deformed ancestor algebra (3.9). 

It is important to note that Lax operators like (3.4) and (3.8) in their bosonic realization appeared 
in some earlier work [10] and they were shown to be the most general possible form in their respective 
class. 

Apart from the discrete models obtained above, one can construct integrable QFT models from 
their exact lattice versions at the continuum limit. For such construction we have to scale the operators 
like Apj, Ac^ , A2^j^ with lattice spacing A, which at continuum limit A — > give pj p{x),ipj — 
tp{x) etc. In this process Lax operator £.{x, A) for the continuum model is obtained from its discrete 
counterpart as Lj{X) — J + AjC{x). The associated i?-matrix however remains the same since it does 
not contain A. Thus integrable quantum field models like sine-Gordon, Liouville, NLS or the derivative 
NLS can be obtained from their discrete variants constructed above. 

5 Quantum integrable inhomogeneous models 

In spite of intensive study of integrable inhomogeneous classical models, their quantum versions seem 
to remain unexplored, except perhaps the recent work [15]. We can construct a novel class of quantum 
integrable models as inhomogeneous versions of the original integrable models listed above. The idea 
of construction is to take the values of the deforming operators to be site dependent functions, which 
would replace in g{u) (3.7) by site dependent operators M^. Consequently all c's in Lax operator 
(3.4) would be site-dependent Cj's. 

However the YB algebra in inhomogeneous lattice models remains the same as in their original 
models with the same quantum i?-matrix. Physical interpretation of such inhomogeneities may be as 
impurities, varying external fields, incommensurate-Ness etc. Few examples of such models are 

1) . Inhomogeneous sine- Gordon model: 

Taking the values of deforming operators different at different lattice sites: cf = Amj-e*"^-', = 
Amje~'"^j , which yields M+ = (Am^)^, M~ = 0, one can construct a novel variable mass discrete SG 
model without spoiling its integrability. In the field limit A — it would yield a quantum integrable 
inhomogeneous sine-Gordon model with variable mass m{x) in an external gauge field 6{x). 

2) . Inhom,ogeneous lattice NLS raodel 

Consider site-dependent (in general also time-dependent) deforming operators in (3.8) and in (4.13) 

as 

ci=(^ = 9j(t), cl = -cl = fj{t), giving m+ = g], m" = 0, (5.16) 

with fj,gj time dependent arbitrary discrete functions. This gives from (4.13) ggiN) = gj{2sj — NjA), 
which reduces (3.8) to the L-operator of the quantum integrable inhomogeneous exact lattice NLS. At 
continuum and high spin limit Sj —>■ -^gj^ limit we get a quantum integrable inhomogeneous NLS field 
model with two arbitrary functions f{x,t),g{x,t) [15]. Note that unlike conventional inhomogeneous 
classical models, our construction is iso-spectral, since inhomogeneity enters here solely through the 
deforming operators in the Lax operator without affecting the spectral parameter and hence the R- 
matrix. 



8 



3). Quantum integrable inhomogeneous Toda chain 

We choose inhomogeneity as c\j,Cij with 03 = 03 = 0, resulting = 0. This gives from (3.8) 
quantum integrable inhomogeneous Toda chain with the Hamiltonian 



In a similar way inhomogeneous versions of Liouville model, relativistic Toda, Ablowitz-Ladik 
model etc. can be constructed. 

6 Quantum integrable hybrid models 

A new class of integrable models may be constructed combining different integrable models which 
interact between themselves and preserve their integrability as a hybrid model. 

The idea of construction is to take Lax operators Lj^^ of different descendant models a = 1, 2, 3.., 
having the same i?-matrix (as listed above) and insert them at different sites j through T(A) = 
L^f^ ■ ■ ■ , which satisfies the QYBE and hence represents a quantum integrable model. Few 
interesting examples are: 

1) . Massive Thirring model : A spectacular example is a new bosonic version of quantum integrable 
massive Thirring model [8, 4]. It can be constructed as a hybrid of two integrable DNLS models by 
fusing two DNLS Lax operators : LoNLSii^^i^i) ® -^B7VLS2(|' ^2) yielding the Lax operator of the 
massive Thirring model involving two bosonic fields (?/'i,V'2)- 

2) . Integrable matter radiation models: Following the above idea we may construct a series of important 
matter-radiation models and their q-deformations as quantum integrable hybrid models [16]. We 
furnish more details in the next section. 

3) . Exotic Hybrid models: We can construct in principle hybrid integrable lattice models by combining 
Lax operators of any lattice models belonging to the same class. For example one can build hybrids 



i) Toda chain and discrete NLS model, ii ) Relativistic Toda chain and the Ablowitz-Ladik model 

iii) lattice sine-Gordon with Liouville model etc. 
4) Integrable quantum field models with defect 

It is challenging to extend the concept of hybrid models to quantum field models. A particular 
direction of this program is to construct integrable quantum field model with defect, e.g., to construct 
NLS, sine-Gordon (SG) or Liouville QFT model with defect. 

Corrigan et al [17] have applied similar idea for constructing SG equation at a; < 0,x > and a 
defect at X = 0. Their approach is concentrated mostly on the soliton scattering problem in analogy 
with the classical case. It is however a challenge to construct and solve a genuine quantum integrable 
SG field model with defect. We could make such a breakthrough very recently by building a quantum 
integrable exact lattice version of the SG model with defect which satisfies the QYBE and allows exact 
solution through algebraic Bethe ansatz [18]. We furnish some details in sect. 8. 




(5.17) 



of 



9 



7 Quantum integrable multi-atom matter-radiation models and their 
q-deformations 

The strategy of construction is to build the monodromy matrix of the system by taking a combination 
of different Lax operators: T(A) = L^{X) Ylf'^ ^jWi with -L*(A) hnked either to the rational (3.8) or 
to the trigonometric (3.4) ancestor model, while A^^-number of Lj{X) are related to the spin or the 
q-spin model. By construction T(A) must satisfy the YB equation with i?''"*(A) or the i?*''*^ (A)-matrix 
in the first or in the second case, respectively, representing a quantum integrable system in both the 
cases with t(A) = tr T{X), generating a commuting set of conserved operators. In our matter-radiation 
models, the subsystem described by L^{X) with the bosonic or the q-bosonic realization represents the 
single mode radiation field, while Lj{X) (in spin or q-spin realization) represents atoms, i.e. the matter 
part of the model. 

Thus through different reductions of a general model we can generate the following integrable 
multi-atom MR models in a unified way [16]. 

1) . Integrable multi-atom Buck-Sukumar (BS) model: 

The model is given hy H = H(i + Hss + Hss with 

Na Na 

Hbs = u;fs' + ^ (uajS] + a{s+Sr + 3-3+)) + a J2iS-S+ - S+S-), (7.18) 

j i<j 

which with a bosonic realization of su(l,l): = \/lVb^,s~ = b\/^,s^ = N + ^ and the spin-s 
operator S = X^^* ak, represents an integrable multi-atom BS model with inter-atomic interactions 
and nondegenerate atomic frequencies (AF). At Na = 1, (7.18) recovers the known BS model [19]. 

2) . integrable multi-atom Jayanes-Cummings (JC) model: 

With bosonic realization s~ = b,s~^ = b^,s^ = b^b, we get a similar multi-atom JC type model. 
The known JC model [20] is recovered at Na = 1, when interatomic couplings vanish and all AF 
coincide. 

3) . integrable trapped ion (TI) model 

with realization = e^*^, = p -|- x, we get a new quantum integrable trapped ion model with 
exponential nonlinearity, having the form (at Na = 1) 

HTi = {ua-Uf)S' + S'^ + a{e-'^S+ + e'''S-) + H^p, H^^ = ]^{p' + x") + xp, (7.19) 
7.1 Integrable q-deformed MR models 

The strategy here is the same, only one has to start now from the trigonometric ancestor Lax operator 
(3.4) and associated i?*^*^-matrix. The Hamiltonian for our q-deformed MR models takes the form 
(for Na = 1) 

HqMR = Hd + {s^Sq + s~Sg ) sina, = —icoCos{aX) + csin(aX), (7.20) 

with X = (s^ — Sg + to), which represents a new class of integrable MR models with € suq(2)) 
and Sq generating (3.5). We list below a new scries of q-deformed integrable MR models which we 
obtain from the same (7.20) for different realizations of Sq. 
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1) . Integrable q-deformed BS model 

This may be constructed from (7.20) at cq = 0, by realizing Sg through q-boson: 

~ \J\.-^\<Mi^ ~ ^q\JV^\q^ = + 5) quantum spin operator Sq through its co-product : 

Sq = Y.jQ *cr|^g^'>^ \ = Yfj^^- Note that at s = 1, we get an integrable version of an 

earlier model [22]. 

2) . Integrable q-deformed JC model 

The model can be constructed similarly from the same general model (7.20) with choice cq = i,c = 1 
and q-bosonic realization s+ = = 6g, = A''. 

3) . Integrable q-deformed TI model: 

Under reduction cq = i,c = and realizing through canonical operators: = e^*^, = p, we 
can construct a q-deformed TI model again from (7.20). 

By taking higher Na values multi-atom integrable variants of the above q-deformed matter- 
radiation models can be constructed. 



8 Quantum integrable sine- Gordon field model with defect 

At the discrete level, as discussed above, it is easy to build hybrid models inserting Lax operators 
of different integrable systems at different lattice sites. However when we try to take the continuum 
limit for constructing the corresponding field model, the above construction breaks down due to the 
lack of any overlapping smooth region. Therefore one needs certain mechanism similar to the classical 
formulation of the same problem [17]. 

We report here about some breakthrough in solving this problem, with details to be presented 
separately [18]. We construct our quantum integrable sine-Gordon model with defect first as an exact 
lattice version with its monodromy matrix defined as 

T\iO = {LU^,y-ti) ■ ■■Lti^,ut)) F^i^,u+,u^) {LZ^i^,uZ,) ■ ■ ■ LZr,{^,uZj,)) (8.21) 

where Lfi^, u^), j = ±1, . . . ± A*" is the quantum Lax operator of exact lattice SG model [5] for the 
field u^. At the defect point j = we define an overlapping quantum Lax operator 

where e± = e^''*?'^"') \ pigt^Po^J'o) , with [uk,Pj\ = i6kj- Remarkably (8.22) is linked with the 
ancestor Lax operator (3.4) under reduction = 0, a = 1, 2, and therefore must satisfy the QYBE 
(1.1) with matrix, which is also the case with the sine-Gordon Lax operators appearing in 
(8.21). This proves the quantum integrability of the system with t(^) = tr{T^j^{^)), generating the 
mutually commuting set of conserved operators. 

It is crucial that this discrete integrable system has a smooth continuum limit: A ^ 0, with the 
fields u^{x),p'^ Ap^{x), having [u^ (x) , p^ (y)] = i6{x — y) and cr^Lf ^ (1 + AU^'^^{x)), 

for x G (0"*", +00) and x G (—00, 0^), respectively, with U^'^^{x) giving the Lax operator for the sine- 
Gordon field model with fields M=''(t, x). At the defect point we get the transition: Fq ^ (l+A£(x))Fo, 
at X e [0~,0+], where C{x) = {Fq^^ + P^){x)Pq^{x): where FP is obtained from F'' by replacing 
P"^ — > |(p~(x) ±p"'"(x)) and Fq gives the known Backliind operator for the sine-Gordon model [23]. 

Starting from the Lax operators at different regions we can find all conserved quantities including 
Hamiltonian, momentum etc. defined on the whole axis, (including the defect point) by using the 
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standard Riccati equation technique. We can also find the corresponding eigenvalues exactly by 
applying the algebraic Bethe ansatz method [24] to the lattice regularised model (8.21) we have 
started with. 

9 Concluding remarks 

We have presented here a unified construction of quantum integrable models from a single ancestor 
model based on a general Yang-Baxter Hopf algebra, representing a new operator deformed quantum 
algebra. Though some of our results were published earlier, we present our whole scheme in a system- 
atic way focusing on a new concept of operator-deformation of an algebra, which is our main concern. 
We also report on some completely new result on the quantum integrable sine-Gordon model with 
defect and other hybrid models like integrable q-deformation of a few matter-radiation models. 

We emphasize that all integrable models listed here allow exact Bethe ansatz solution. Moreover 
, similar to their unified construction one can get their solution also in a unified and almost model- 
independent way. For this one has to start from the general sohition of the ancestor model following 
the algebraic Bethe ansatz and subsequently obtain the results for all descendant models with the 
same i?-matrix as various reductions of the general result. 
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